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ABSTRACT 

We  study  the  system  of  equations  (1.1),  describing  the  stationary  motion 
of  a  compressible  viscous  fluid  in  a  bounded  domain  ft  of  R3.  The  total 
mass  of  fluid  mjft|,  inside  ft,  is  fixed  (condition  (1.2)).  We  prove  that  for 
small  f  and  g,  there  exists  a  unique  solution  (u,p)  of  the  above  system 

of  equations,  in  a  neighborhood  of  (0,  m).  Moreover,  by  introducing  a 

(«■-'  (c  -  ^ 

suitable  parameter  X,  we  prove  that  the  solution  of  the  Navier-Stokes 
equations  (<U~t49  are  the  incompressible  limit  of  the  solutions  of  the 
conpressible  Navier-Stokes  equations  .  The  proofs  given  here,  apply, 

without  supplementary  difficulties,  in  the  context  of  Sobolev  spaces  H^'J^, 
and  other  functional  spaces.  The  results  can  be  extended  to  the  heat 
depending  case,  too. 

AMS(MOS)  Subject  Classifications:  35G30,  35M05,  35Q10,  76D05,  76N10 
*  Key  words:  Non-linear  partial  differential  equations,'  Viscous 


compressible  fluids;  Incompressible  limit.  Stationary  solutions 


Work  Unit  Number  1  -  Applied  Analysis 


/* 


^Department  of  Mathematics,  University  of  Trento  (Italy). 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


SIGNIFICANCE  AND  EXPLANATION 


c 

In  this  paper  we  considers  the  non-linear  system  of  partial  differential 

equations  (1*lf,  describing  the  barotropic  stationary  motion  of  a  compressible 

,  •  '  " k  Y' 

fluid,  in  a  bounded  region  Os  see  [4]  for  details. 

We "assume  that  the  total  mass  of  fluid  inside  Q  is  fixed,  and  equal  to 
mjn|,  where  the  mean  density  m  is  given,  (condition  (1.2)). 

We  prove  That  ^or  small  f  and  g,  there  exists  a  unique  solution  u(x), 
p(x)  (Llf' in  a  neighborhood  of  (0,  m).  Here,  u(x)  is  the  field  of 

velocities,  p(x)  the  density  of  the  fluid,  p(p(x))  the  pressure  field, 

'  f  >  W  t  *  c 

and  f(x)  the  external  force  field  (in  the  physical  interesting  case  one 
has  g  *  0). 


Moreover,  wo  prove  th*<j  the  solutions  of  system  converge  to  the 

solution  of  the  Navier-Stokes  equation  ,'  as  X  +•;  i.e.  when  the  Mach 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MEG,  and  not  with  the  author  of  this  report. 


STATIONARY  MOTIONS  AND  INCOMPRESSIBLE  LIMIT  FOR  COMPRESSIBLE  VISCOUS  FLUIDS 

H.  Beirao  da  Veiga* 


1.  INTRODUCTION  AND  MAIN  RESULTS 

In  this  papsr  we  study  the  systsm 

-yAu  -  vVdiv  u  +  7p(p)  ■>  p[f  -  (u»V)u] ,  in  ft, 

(1  div(pu)  -  g,  in  ft, 

[  u|r  "  °' 

in  a  bounded,  open  domain  in  R3,  locally  situated  on  one  side  of  its  boundary  T,  a  C3 
manifold  (the  case  n  |>  3  can  be  studied  by  the  same  method).  As  usual, 


<*•*>«-  l  vi  IT  ■ 

i-1  1  8 l 


System  (1.1)  describes  the  stationary  motion  of  a  barotropic,  compressible  fluid;  see 
Serrin  (4).  In  equation  (1.1),  p(x)  is  the  density  of  the  fluid,  u(x)  the  velocity 
field,  f(x)  the  assigned  external  force  field,  p  -  p(p)  the  pressure.  In  the  physical 
equation  one  has  g  *  0;  however,  on  studying  (1.1)  from  a  mathematical  point  of  view,  it 
is  not  without  interest  to  study  the  general  case. 

We  assume  that  the  total  mass  of  fluid  inside  0  is  fixed,  i.e.,  we  impose  to  the 


solution  of  (1.1)  the  contraint 


w  fn  oU,dx  “ 


where  the  mean  density  m  is  a  given  positive  constant.  The  function  p  will  be  written 
in  the  form  p  »  a  +  a,  and  the  new  unknown  o(x)  has  to  verify  the  constraint 


v  w, 
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.  -  .’-Ai 


(1.3) 


5  5  w  a(x)dx  ■ 0 


We  assume  Chat  the  real  function  p  *  p(p)  is  defined  and  has  a  Lipschitz  continuous 
first  derivative  p'(p)  In  a  neighborhood  I  =  [m  -  l,  m  +  i)  of  m,  for  some  positive 
l  <  m/2.  We  assume  also  the  (unessential)  physical  condition  k  -  p* (m)  >  0.  Clearly, 
(1.4)  p'(p)  “  p' (m  +  c)  «  k  -  u(o),  V  o  €  I, 

where  u>(c)  is  a  Lipschitz  continuous  function,  such  that  w(0)  «  0.  We  set 


„  _  .  lute)  -  «i(t)  I 

S  —  BUD  1  (  ■  '  11  I  “  e 

a  -  t 

a, tel  1  ' 

Concerning  the  constants  u  and  v,  we  only  assume  that 
(1.5)  u  >  0  ,  v  >  -u. 

In  the  sequel,  we  write  the  system  (1.1)  in  the  equivalent  form 

/ 

-pAu  -  vVdiv  u  +  kVa  «  w(a)  Va  + 

+  (m  +  a)  [f  -  (u*V)u),  in  £), 
m  div  u  +  u  •  Vo  +  a  div  u  ”  g,  in  (2, 

U|r  -  °  • 

Let  us  introduce  some  notation.  We  set 


Vv 


3  3v 

I 

i,k=1  k 


Vv 


V2t 


3  3v 

I  3 r 

i,k=l  k 


32t 

3V*k 


where  v  is  a  vector  and  T  a  scalar. 

We  denote  by  Hk,  k  integer,  the  Sobolev  space  W k,2(Q),  endowed  with  the  usual  norm 

I  1^,  and  by  I  |p,  1  <  p  <  +“,  the  usual  norm  in  L^  «  L^(0).  Hence,  I  1q  “  I  Ij*  For 

convenience,  we  utilize  the  same  symbol  H*1  to  denote  also  the  space  of  vector  fields 

k  2 

v  in  ft  such  that  r  W  '  (W#  i  ■  1#  2,  3.  This  convention  applies  to  all  the 
functional  spaces  and  norms  utilized  here. 


V-.-VV  Vk-w  ■/«  '.-V.  VT  WlV.VA'l'.'A’Vr^, '  M  v.  w*.  ^  *.■  *.*.>  A  i_M_< ,  •  1 


Ti’V"."’ 


For  k  >  1,  we  define 

h£  “  {v  e  H* :  v  »  0  on  T)  . 

Moreover , 

H2  -  {t  €  H2s  t  -  0}  ,  H2  -  H2  D  H2, 

where  t  is  the  nean  value  in  fl  of  the  scalar  field  t(x).  Finally,  for  vector  fields, 
we  define 

H2  .  «  {v  t  H2  :  div  v  «  0  on  T}. 
u ,  a  u 


In  the  sequel,  c,  cQ,  c.,,  c2, • • • ,  denote  positive  constants  depending  at  most  on  0. 
Moreover,  c',  eg,  c^,  . ..,  denote  positive  constants  depending  at  moat  on 
fl,  y,  v,  k,  m,  t,  and  S.  The  same  symbol  c  (or  c* )  will  be  utilised  to  denote 
different  constants,  even  in  the  same  equation. 

In  section  3  we  prove  the  following  result: 

Theorem  *.  There  exists  positive  constants  eg  and  c  ^  such  that  if 
1  -2 

f  e  H  ,  g  f  Hfl,  and 

(1.7)  Ifl,  ♦  lgl2  <  cQ. 

3  -2 

then  there  exists  a  unique  solution  (u ,o)  e  HQ  x  H  of  problem  (1.6),  in  the  ball 


(1.8) 


lul3  +  lolj  <  Cj. 


K  crucial  tool  in  order  to  prove  this  result  will  be  the  study  of  the  linear  system 

-yAu  -  vVdiv  u  +  kVo  •  F,  in  S), 

)  m  div  u  +  v  •  Vo  +  o  div  v  ■  g,  in  {5, 

u|r  "  0  • 

for  which  we  will  prove  the  following  result: 


-3- 


3  -2 

holds.  Then,  there  exists  a  unique  solution  (u,c)  t  HQ  ^  *  K  of  the  linsar  svst— 
(1.9).  Moreover, 


(1.10) 


plul,  +  klol,  <  c(l  +  >*  * +  -li-i- 
3  2  v  g  +  v  s 


In  section  4,  we  assume  that  the  function  p(p,X)  depends,  in  a  suitable  way,  on  a 
parameter  1.  By  letting  X  *  +»,  we  prove  that  the  solution  of  the  Navier-Stokes 
equation  (1.14)  is  the  incompressible  limit  of  the  solutions  of  system  (1.13).  For  the 
justification  of  the  physical  aspects  of  the  description  (i.e.,  the  behavior  of  p(p,X), 
as  X  ♦  +<■>)  we  refer,  for  instance,  to  reference  [1]. 

We  assume  that  for  each  value  of  the  parameter  X  e  [XQ,  +  <•[  (XQ  e  It,  has  no  special 
meaning)  the  function  p(p,X)  is  defined  in  a  neighborhood  £  [m-l^,  m+t^l  of  m, 
where  0  <  tx  <  m/2.  Moreover,  for  each  fixed  X,  the  derivative  dp(p,X)/dp  £  p'  (p,X), 
is  Lipschitz  continuous  on  1^,  with  Lipschitz  constant  S^. 

We  define  k^  £  p'(m,X),  and  assusw  tht  k^  >  kQ  >  0  (the  constant  kQ  has  no 
special  meaning,  since  we  will  let  k^  ♦  -Hi,  as  X  ♦  +*•) .  We  suppose  that  there  exist 
positive  constants  4  and  1  such  that 

(1.11)  <  ♦  kf  ,  V  X  >  X„  , 
and 

(1.12)  ixk^  >  t,  V  X  >  XQ  . 

By  (eventually)  defining  a  smaller  t^,  we  assume,  without  losing  generality,  that 
i^k^  -  t.  Finally,  let  <1)^(0)  be  defined  by  p'(m+o,X)  »  k  -  u^to). 


Consider  the  stationary  compressible  Navier-Stokes  equation,  with  state  function 


.  1  Mil. 


m  v  ■  71  p*7  >,J  y.'yr 


31!.  i*. 


w.v-r*: 


^VPC  iNi-f 


I  • 


I 


i 


(1.13) 


-U4ux  -  v7  dlv  +  7p(p^,X)  ■  p^[f  -  (u^.VJu^l  , 
divlp^u^)  “0  ,  in  ft  , 


{ux> 1 r  “  0 


and  the  incompressible  Navier-Stokes  equation 

~y4u^  +  7ir(x)  “  m[f  -  , 

(  div  us  ■  0  ,  in  ft  , 


(1.14) 


<UJ| 


As  above,  we  set  p^(x)  “  m  +  o^(x),  and  we  look  for  solutions  of  (1.13)  verifying 
assumption  (1.2),  i.e.  such  that  (1.3)  holds. 

He  denote  by  c,  cQ,  ,  c.,...,  positive  constants  depending  at  most  on 

ft,  u,  v,  m,  1,  *  and  kQ,  and  we  say  that  a  positive  constant  is  of  type  c  if  it 
depends  at  most  on  the  above  parameters. 

In  section  4,  we  prove  the  following  result i 

Theorem  C.  There  exists  positive  constants  cQ  and  such  that  the  following 

statement  holds-. 

(i)  Let  f  f  H1,  belong  to  the  ball 
(1.15)  If  I.  < 


'1  *  -0  ’ 


3  _  -2 


Then,  for  each  X  >  X0,  problem  (1.13)  has  a  unique  solution  (u^,o^)  t  HQ  x  H  in 


the  ball 


(1.16) 


* UX ' 3  +  VV2  *  C1 


(ii)  If  lim  k.  «  +«,  then 


-5- 


m. 


.  •  .  '  *  1 


\  \  V 


-?l  ^  iT-  n. 


(1.17) 


♦  Mm  ,  weakly  in  Hq,  strongly  in  H®,  V  s  <  3  , 
div  u.  ♦  0  ,  weakly  in  H«  ,  strongly  in  H5,  W  s  <  2  , 


♦  0  ,  strongly  in  H  , 

Vp(p^,X)  +  Vir  ,  weakly  in  H1,  strongly  in  H®,  V  s  <  1  , 

where  (u^.Vx)  is  the  unique  solution  of  problem  (1.14). 

The  existence  of  the  solution  (u^.ViO  of  (1.14)  is  well  known.  However,  it  follows 
from  our  proof,  too. 

Note  that  both  problems  (1.13),  (1.14)  are  invariant  under  addition  of  arbitrary 
constants  to  p(X,p)  and  n,  respectively. 

An  existence  result  for  system  (1.1)  was  given  first  by  Padula,  in  reference  (3). 
Unfortunately,  the  (quite  simple)  proof  given  there  depends  in  a  crucial  way  on  a 
smallness  condition  on  p  respect  to  v  (  g  and  v  positive  constants).  This 
condition  was  dropped  in  Valli's  paper  (5),  where  a  result  similar  to  Theorem  A  is  proved, 
by  approximating  the  stationary  solutions  with  the  periodic  solutions  of  the  corresponding 
evolution  problem.  This  technique  was  applied  in  [6]  to  the  heat-depending  case,  and  to 
more  general  boundary  conditions. 

The  proofs  given  in  our  paper  are  quite  simple,  and  apply  as  well  (without  any 
supplementary  difficulty)  in  the  context  of  other  spaces  of  functions,  as  for  instance 
Sobolev  spaces  H*5'*3,  1  <  p  <  +«>. 

In  particular,  for  small  data  (f,g)  e  Hk+1  x  2 ,  there  exists  a  unique  solution 
k+3  — k+2 

(u,o)  t  Hq  x  h  ,  in  a  neighborhood  of  the  origin,  for  every  k  >  0  (we  assume  the 
derivative  p*k+1^  Lipschitz  continuous,  and  0  of  class 

Furthermore,  all  the  results  hold  again  in  any  dimension  of  space  (on  dealing  with 
the  non-linear  problem  in  space  H^'P,  k  must  be  sufficiently  large). 

Statements  and  proofs,  in  the  above  general  setting  up,  will  be  given  in  a 
forthcoming  paper,  where  (for  completeness)  we  will  consider  the  heat-conductive-case.  In 
this  paper,  we  state  only  the  counterparts  of  theorems  A  and  B,  in  appendix  2.  The  proofs 


& 


can  b«  easily  done,  by  following  those  of  theorems  A  and  B.  Here,  we  have  preferred 


considering  the  main  case  (1.1)  by  itself,  in  order  to  avoid  secondary  technicalities.  In 
fact,  in  the  heat-depending  case  a  third  equation  should  be  added  to  system  (1.1)  (see 

(6.1) ).  However,  that  equation  is  weakly  coupled  with  its  companion  equations,  in  system 

(6.1) .  As  a  matter  of  fact,  the  more  interesting  mathematical  problesu  and  the  main 
difficulties,  already  appear  on  studying  system  (1.1). 

Finally,  it  goes  by  itself  that  quite  obvious  assumptions  and  devices,  allow  the 
coefficients  u»  v  x#  «v<  "V  «nd  y*  on  depending  on  u,  p  and  6. 

2.  Proof  of  Theorem  B 

We  start  by  proving  the  uniqueness  of  the  solution  of  the  linear  system  (1.9),  under 
the  assumption  (2.1)  below.  Let  (u,o)  be  a  solution,  with  data  F  «  0,  g  »  0.  By 
multiplying  both  sides  of  equation  (1.9)^  by  mu  and  of  equation  (1.9)2  by  k a,  by 
integrating  over  fl.  and  by  adding  side  by  side  the  two  equations,  one  easily  shows  that 

muQ  «Vulg  <  |  |div  vl,  loljj  , 

where  UQ  ”  min{y,  u  +  v}. 

Hence , 

lul^  <  c  — —  Idiv  v|  lali!  . 

1  m^  1  1 "  0 

Moreover,  from  (1.9)  ^  it  follows  that 

klolg  <  cklVol_.j  <  C(u  ♦  |v|)  lul1# 

since  a  “  0.  Consequently, 


<U  »  M> 


|div  vj^  lul,. 


-7- 


for  a  suitable  positive  constant  cqi  recall  that  H  C  L 


In  the  remaining  of  this  section  we  prove  the  existence  of  the  solution  of  system 


(1.9).  He  assume  that  v  e  H„  .  verifies  the  condition 

0,  <3 


(2.2) 


|div  v|M  +  2j Vvl^  <  rak/(y  +  v)  , 


1  -2  -2 

and  that  F  r  H  ,  g  e  H  •  let  t  e  H  ,  and  consider  the  linear  problem 


(2.3) 

where 


mk 


y  +  v 


X  +  v  •  VX  “  G  , 


div  F  - 


(2.4) 


G  «  Ag  +  — 7 — 

’  y  +■  v 

-  [2Vv  :  V^T  +  Av  •  Vt  +  A  (r  div  v)]  . 


The  significance  of  equation  (2.3)  is  strongly  related  to  the  identity  (2.20).  It  is 

2 


well  known  (Lax-Phillips  (2) )  that  there  exists  a  linear  map  G  ♦  X,  from  all  of  L' 

5  2 

into  L  ,  such  that  for  each  G  e  L  the  corresponding  X  is  a  weak  solution  of  (2.3), 
and  verifies  the  estimate 
(2.5) 


4  —  IXI„  <  1Gln. 

2  u  +  V  0  0 


By  a  weak  solution  of  (2.3),  we  mean  here  a  function  X  e  L  such  that 


(2.6) 


— — —  f  X»dx  -  f  X  div(»v)dx  =  /  G»dx,  V<p  e  H1 

u  +  v  n  1  n  n 


For  the  reader's  convenience,  >»«nplete  proof  of  this  result  in  the  Appendix 

I. 

14  2  <» 

By  using  the  embeddings  H  CL  and  H  CL,  one  verifies  that  IGIQ  is  bounded 
by  the  right  hand  side  of  equation  (2.7)  below.  Hence  our  solution  X  of  (2.3)  verifies 


(2.7) 


g  +  v 


IXIQ  <  cllF^  +  I  g  1 2  +  lvl3  ItIj) 


Let  now  9  f  H  be  the  solution  of  the  Dirichlet  problem 


(2.8) 


(g  +  v)  A9  =  kX  -  div  F,  in  ft 


e|r  '  0  • 


By  using  (2.7),  one  has 


(2.9) 


(g  +  v) 


161  <  C  —  ((FI  +  1  g  1  +  Ivl  I T  1  _  )  +  C I  FI  . 

2  m  l  2  5  2  i 


Define  now 

(2.10) 


eQ(x)  =  9(x)  -  9  . 


Clearly,  9Q  =  0.  Let  (u,  a)  be  the  unique  solution  in  the  class 
following  linear  Stokes  problem,  in  S3  : 


(2.11) 


S 

-gAu  +  k?o  =  F  +  vV6p, 
/  div  u  «  6q, 


From  the 


estimates  for  this  problem  one  has 


(2.12)  g«ul3  +  klol2  <  cdFl,  +  |  v  |  19  «2  +  glB^)  . 

By  taking  in  account  that  *  *®!2'  one  9e*-s 

glul  +  klal  <  cfl  +  jf«  + 

3  2  y  +  v  m  l 

(2.13) 

.NS 

1  m 


+  c. 


( Igl,  +  ivl  +  Irl,)  . 


IX 


Let  now  c2  be  a  positive  constant  such  that 

! div  w|„  +  2j yw)^  <  c2lwl3,  Vw  e  H2  . 

In  the  remaining  of  this  section  we  assume  that  the  vector  field  v  verifies  the 
condition 


(2.14) 

where,  by  definition, 
(2. 14') 


lvl3  <  yk  , 


y  s  mini 


u0» 


c0(m>M)2  '  2c,<Hv|)  *  c2<u+v)> 


Assumption  (2.14),  implies,  in  particular,  (2.1)  and  (2.2). 

From  (2.13)  and  (2.14),  one  gets 

uJuJ3  +  klo«2  <  k  It  1 2  + 

(2.15) 

♦  c(l  ♦  JL+-M  ♦  ,f,  +  c  jl±JvL  ,gl  . 

'  U  +  v  m  '  1  m  3  2 

At  this  point,  we  call  attention  to  the  sequence  of  linear  maps,  introduced  above: 

(F,  g,  T)  +  (F,  \)  *  (F,  0)  +  (F,  0Q )  ♦  (u,  a)  , 
which  were  defined  by  equations  (2.3)  +  (2.4),  (2.8),  (2.9),  (2,11),  respectively.  The 
product  map  (F,  g,  t)  *  (u,  a)  is  linear  and  continuous,  by  (2.15).  Hence,  if 
(u1#  o^)  is  the  solution  corresponding  to  data  (F,  g,  t),  it  follows  that 
(u  -  Uy  o  -  ct ^ )  is  the  solution  corresponding  to  data  (0,  0,  t  -  T ^ )  .  Consequently, 
(2.15)  yields,  in  particular, 

la  -  o,l2  <  }  «t  -  Ttl2  . 

-  2 

Hence,  for  fixed  F  and  g,  the  map  t  +  a  is  a  contraction  in  H  .  Consequently 
it  has  a  (unique)  fixed  point  a  ”  T  . 

In  the  sequel  we  prove  that  the  pair  (u,  o),  corresponding  to  the  fixed  point 
o  ”  t,  solves  equation  (1.9).  Equations  ( 1 . 9 )  |  and  (1.9)j  follows  from  (2.11).  In  order 


to  prove  ( 1 . 9 ) 2 ,  we  start  by  substituting  the  expression  of  obtained  from  equation 
(2.8).,  in  the  first  term  on  the  left  hand  side  of  (2.3).  This  yields,  since  t  “  a  , 


(2.16) 


mA9  +  v  •  7X  +  27v  t  7  o  *  Av 


7o  +  A (c  div  v)  » 
«  Ag  • 


On  the  other  hand,  by  applying  the  divergence  operator  to  both  sides  of  equation 
(2.11),,  and  by  utilizing  (2.11)2  one  gets  -{y  +  v)A0  +  kAo  “  div  F,  since  A90  -  A8.  By 
comparison  with  (2.8),,  one  shows  that  X  “  Ac.  By  replacing  X  by  Ao  in  equation 

(2.16) ,  it  follows  that 

2 

mA  div  u  ♦  v  •  7A c  ♦  27v  i  7  c  +  Av  •  7c  + 

(2.17) 

♦  A(c  div  v)  -  g  »  0  , 


or  equivalently. 


(2.18)  A[m  div  u  +  v  •  7o  +  c  div  v  -  g)  ■  0,  in  ft  . 

The  function  between  square  brackets  (which  belongs  to  H1)  is  equal  to  the  constant 
-m9  on  the  boundary,  by  (2.8)2,  (2.10),  (2.11)2,  and  by  the  assumptions  v  •  0, 
div  v  »  g  =>  0  on  T  .  Consequently, 

m  div  u  v  •  7o  +  c  div  v  -  g  »  -ra9,  in  ft  . 

By  integrating  both  sides  of  this  equation  in  ft,  one  shows  that  it  must  be  0  «  0  . 
Hence,  equation  (1.9)2  is  satisfied.  Finally,  the  estimate  (1.10)  follows  from  (2.15). 

a 

Remark.  One  has  to  be  careful  on  deducing  (2.18)  from  (2.17),  since  both  equations  hold 
only  in  a  weak  sense.  The  point  is  to  prove  the  identity 

“/ft  *  ?fdx  “  -/  Ac  div(vv)dx  + 


(2.19) 


/„  [27v 


72o 


♦  Av 


7o]  edx,  V*  t  cj  , 


which  is  a  weak  formulation  of 

(2.20)  A(v  •  7c)  “  v  •  7Ao  +  27v  :  72o  +  Av  •  7o  • 
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For  a  €  H  ,  this  last  identity  holds,  and  yields  (2.19)*  If  o  t  H  ,  we  approximate 

2  3 

it  (in  the  H  norm)  by  a  sequence  of  functions  on  €  H  ,  and  we  pass  to  the  limit  in 

equation  (2.19)  (written  with  o  replaced  by  a  )  as  n  ♦  ♦  •  • 

n 

3.  Proof  of  Theorem  A 

For  convenience,  in  this  section  we  will  not  take  case  on  the  explicit  dependence  of 
the  positive  constants  respect  to  the  parameters.  Bowever,  all  the  constants  depend  at 
most  on  ft,  u,  v,  k,  a,  l  and  S. 

_2 

Let  c3  be  a  constant  such  that  | T | ^  <  c3  ItIj,  for  every  t  «  H  .  we  will 
utilize  here  the  condition 


which  guarantees  that  m  +  a(x)  belongs  to  the  domain  of  p,  for  every  x  e  fi,  since 
-l  4  t(x)  4  t. 

Let  v  (  h|  verify  (2.14),  and  T  e  H2  verify  (3.1),  define 
(3.2)  F( v,  T)  “  (m  +  x)(f  -  (v  •  V)v]  “  u(t)Vt, 

and  consider  the  linearized  system  (1.9)  with  F(x)  given  by  F(v,  t),  l.e.,  the  system 


(  -yAu  -  vV  div  u  +  kVa  «  F(v,  t),  in  ft, 

<  m  div  u+(v»V)a  +  a  div  v  «  g,  in  (1, 


Since  H1  C  L4,  H2  C  L  and  |w(t>|„  4  S |  x |  ^  <  (t/c^)  S  I  x  '  one  easily  shows  that 
(3.4)  IF(v,  T)l,  4  c<4  m  +  •£-)  (|f|,+  |v»2)  +  cSItI2  . 

1  2  Cj  1  2  2 

This  last  estimate,  together  with  (1.10),  yields  the  following  result: 

Theorem  3.1.  Let  v  e  d,  T  «  H2  and  let  (2.1),  (2.14)  and  (3.1)  be  satisfied. 
Then,  the  unique  solution  u,  a  e  H2  x  H2  of  system  (3.3),  verifies  the  estimate 

lul  +  lot  4  a(  ItI  +  Ivl  )2  +  b(lfl  +  Igl  )  , 


(3.5) 


-.-m  oarj  rji  rji  rm’^n 


where  the  positive  constants  a  and  b  depend  only  on  £5,  y  >  v,  k,  a,  i  and  S. 

The  existence  and  uniqueness  of  the  solution  (u,  o)  of  system  (3.3),  enables  us  to 
define  the  corresponding  map  (u,  a)  ”  T(v,T).  The  fixed  points  of  the  map  T  are  just 
the  solutions  of  the  non-linear  system  (1.6).  In  order  to  prove  the  existence  of  these 
fixed  points  we  assume  that 


(3.6) 
and  that 

(3.7) 


•f, +  '*'2  <-kmin{lr  i;}' 

Ivlj  +  I  t  1 2  <  min{-|^,  yk,  •£-}  . 


The  parameter  6  e  ]0,  1],  will  be  fixed  later  on.  Consider  the  ball 

Bj  I  {(v,  t)  <  H2  ^  x  H2:  (3.7)  holds}  . 

1  2 

This  is  a  compact  set  in  x  L  .  Moreover,  by  using  (3.S),  one  shows  that 

TB.  CB„  for  every  6  <  1  .  We  want  to  prove  that,  for  a  sufficiently  small 
6  6 

6,  depending  only  on  £),  y,  v,  k,  m,  t  and  S,  the  map  T  is  a  contraction  in  Bj . 

Hence  T  has  a  (unique)  fixed  point  in  B ^ ,  and  Theorem  A  is  proved. 

Let  (u,  a)  -  T( v,  T),  (u  1 ,  o1)  «  T(v1#  t ^ )  ,  F  “  F(v,  x),  F1  -  Ffv^  x^  .  One  has, 
in  n, 

f 

-y  i(u  -  u1 )  -  vV  div(u  -  )  +  k7(o  -  o1>  -  F  -  F ^ , 

e)  <  m  div(u  -  u.,)  +  v  •  7(o  -  cr., )  +  (v  -  v1 )  •  7a.,  + 

+  o1  div(v  -  v, )  +  (o  -  a,)  div  v  =  0  . 


By  multiplying  both  sides  of  equation  (3.8),  by  m(u  -  u,)  and  both  sides  of 
equation  (3.8)2  by  k(o  -  c,),  by  integrating  in  £5,  and  by  adding  side  by  side  the  two 
equations  obtained  in  that  way,  one  shows  that 


-13- 


Similarly 


and 


l(v  •  7)v  -  (v1  •  V)v1l_1  <  c(lvl1  +  lv1l1)  Iv  -  v1ll|  , 

It(v  •  7)v  -  t1(v1  •  <  cltlj  I ( v  •  7)v  -  (v1  •  7)v1l_1 

♦  clv,^2  It  -  t 


+ 


ro  • 


Furthermore , 


t(x) 

Iw(t)  7t  -  uIt,)  Vt1«_1  -  17  /  u(C)dCI  < 

X  T  ^  (  X) 

t(x) 

<  1/  *<5>dCl0  <  sdtl.,  ♦  |t, |„)  It  -  T,l  . 

T,(x) 


The  above  inequalities  yield  (3.13). 

For  X  <  a/c^cj,  one  has  c^jdiv  v|b  <  o^e1  Ivlj  <  -j  ,  by  (3.7).  Hence,  from  (3.12), 
(3.13)  one  gets 


(3.14) 


lu  -  u,l2  ♦  lo  -  0,.2  <  c'lfl2  It  -  T  ^  I  q  + 

+  c'[0  +  ItIj)2  ( Ivl ,  +  Iv,!,)2  +  la^l  -  v.,12  + 

+  c[lv,l2  +  S(ItI2  +  It^j)]2  IT  -  t ^ 1 2  . 


By  choosing  S  sufficiently  small,  depending  only  on  0,  v,  v,  k,  m,  l  and  S,  one 


has 


lu  -  Ujl 1  +  lo 


Vo 


<  t(,v 


Vi  +  ,T  ■  Vo) 


Hence  T  is  a  contraction  in  B»,  which  proves  Theorem  X. 


1  2 

Remark.  B ^  ia  a  compact  and  convex  aubaat  of  Bj  «  l  ,  T  i  la  oontlnuoua 

raapact  to  that  topology,  and  TB1  C  .  Hanca,  we  can  prove  the  exlatenoe  of  (at  leaat) 
a  fixed  point  in  Bj  by  using  Shauder'a  theorem.  The  uniqueneaa  follows  by  using  (3.14) 
(actually,  it  ia  quite  trivial  to  obtain  more  atringent  uniqueneaa  reaulta.). 

4.  Proof  of  Theorem  C.  During  the  proof  of  part  (i)  of  theorem  C,  o^,  k^,  «x»  will 

be  denoted  by  u,  o,  k,  u  reapectively.  Theorem  B  atatea  that  if  P  t  H1,  g  »  0,  and  if 
v  c  Hg(<J  verifies  the  condition 

(4.1)  lvl3  <  yk  , 

then  there  exists  a  unique  solution  (u,o)  e  Hg^  *  H2  of  the  linear  system  (1.9). 
Moreover, 

(4.2)  uiuij  kioi2  *  c(i  ♦  +  islbi)  in,  . 

—2 

Let  us  now  fix  t  e  H  in  the  ball 

i 

(4.3)  lrl_  <  —  ,  or  equivalently,  IkTl,  <  —  , 

2  C3  2  C3 

where  c3  was  defined  in  section  3,  and  1  is  the  positive  constant  defined  in  (1.12). 

Condition  (4.3)  guarantees  that  |t(x)|  <  i^,  V  x  e  ft.  In  particular  m/2  <  m  ♦  t(x)  < 

<3m)/2. 

By  defining  F(v,t)  as  in  (3.2)  (recall  that,  now,  u  -  u^)  one  haa,  as  in  section 
3, 

IP(V,T)  I1  <  c(|  m  +  -^(ifl,  ♦  Ivl2)  +  csxlxl2  . 

Hence, 

(4.4)  IP(V,T)1,  <  c(j  a  +  -A_)(|fi1  ♦  Ivl2)  +  c*k2lTlj  , 

(recall  k  -  k^).  If  v  and  t  verify  assusgitions  (4.1)  and  (4.3),  it  follows  from  (4.2) 
and  (4.4)  that  the  unique  solution  (u,o)  of  system  (3.3)  verifies  the  estimate 

(4.5)  lulj  +  Ikol 2  «  a(lkTl2  +  Ivlj)2  +  blf«1  , 

where  now  a  and  b  are  constants  of  type  c  (the  above  result  corresponds  to  theorem 


3.1,  in  section  3). 


The  proof  goas  on  as  in  (action  3,  by  utilicing  now  kc  and  kT  instead  of  o  and 
x,  respectively  (in  this  way,  inequalities  (3.S)  and  (3.1)  becomes  (4.5)  and  (4.3)2, 
respectively t  condition  (4.1)  remains  unchanged). 

Following  section  3,  we  denote  by  T  the  map  (u, o)  ■  T(v,t),  where  the  data 
( v, x )  e  Hgd  x  H2  verify  (4.1),  (4.3),  and  (u,a)  is  the  (corresponding)  solution  of 
system  (3.3). 

He  fix  f  e  H1  verifying  (3.6)  (here,  g  •  0),  and  we  consider  the  restriction  of 


T  to  the  ball  Bj,  0  <  6  £  1,  defined  by  the  condition 
(4.6)  »vl3  «■  Ikrlj  <  minf|j-,  • 


The  substitution  of  t  by  kx  transforms  (3.7)  on  (4.6).  Arguing  as  in  section  3, 


and  recalling  that  k  >  k0,  we  prove  inequalities  (3.10),  (3.11)  and  (3.12),  if  in  these 
inequalities  we  replace  c,  a x,  1 1  by  kc,  kcr^,  kr,  kt^  respectively.  The  constants 
c',  cj,  c‘3,  are  now  of  type  c,  hence  independent  of  k. 

Inequality  (3.13)  holds,  as  written  in  section  3.  Recalling  that  <  4k2,  and 
that  k  >  kg,  we  show  that  (3.13)  holds  again,  if  T,  x},  and  S  are  replaced  by  kx, 
kT.,,  and  *,  respectively,  and  if  the  right  hand  side  of  the  inequality  is  multiplied  by 

1  +  (1/kg). 

By  choosing  5  as  in  section  3,  i.e.  6  <  a/tc^c^),  we  get  an  inequality  similar  to 
(3.14),  where  now  x,  a,  x,  and  0,  are  multiplied  by  k,  and  the  constants  are  of  type 
c.  By  choosing  6  sufficiently  small  (depending  only  on  0,  V »  v,  m,  £,  $,  kg)  one  gets 

(4.7)  ,u"u1,1  +  ,kt’"lt01,g  <  -J  dv-Vjl2  ♦  Ikx-kx^l 2)  . 

Hence  T  is  a  contraction  in  Bj,  which  proves  the  first  part  of  theorem  C. 

We  now  prove  part  (ii)  of  that  theorem.  Condition  (1.15)  guarantees  the  uniqueness 


of  the  solution  of  problem  (1.17),  for  a  sufficiently  small  Cg. 
Let  us  write  system  (1.13)  in  the  form  (1.6),  i.e. 


-17- 


-pAu^  -  vV  div  +  k^Va^  “  “x^°X^aX  +  *  (u^*7)u^]  ' 

m  div  u^  *  u^  •  Va^  +  a ^  div  u^  -  0,  in  0  , 

(u.) I.  =  0  . 


From  (1.16),  it  follows  that  there  exists  u„  e  Hr-',  such  that  (1.17)]  holds  (here, 
we  consider  subsequences  of  u^i  the  convergence  of  all  the  u^  to  u„ ,  as  X  ♦  +•* 
will  follow  from  the  uniqueness  of  the  limit  u„ ,  since  we  will  show  that  u„  is  the 
solution  of  (1.14)). 

The  bound  (1.16),  and  the  hypothesis  k^  ■*  +“  as  X  +  +*,  imply  (1.17) j. 

Furthermore,  equation  (4.8>2,  together  with  (1.16)  and  ( 1 . 17 ) j,  shows  that  div  u^  0, 
strongly  in  Hq,  as  X  +  +”.  Since  Idiv  u^«2  is  bounded,  (1.17)j  follows.  In 
particular,  div  u„  =  0. 

Now,  we  pass  to  the  limit  in  equation  (1.13)],  as  X  +  +“•  One  has  yAu^  ♦  pAu^  and 
-vV  div  ux  +  0,  weakly  in  H1  and  strongly  in  H8,  0  <  s  <  1j  and  px  ♦  m,  strongly 

in  H2.  Moreover,  p^(u^-7)u^  *  “(u^'Vlu^,  weakly  in  H2  and  strongly  in  H8,  0  <  a  <  2. 

By  using  equation  (1.13)],  it  follows  that  Vp(p^,X)  ♦  pAu^  +  mlf  -  (u^VJu^) ,  weakly 

in  H1,  strongly  in  Hs,  0  <  s  <  1.  Obviously,  the  limit  function  must  be  of  the  form 

Vn(x).  Theorem  C  is  completely  proved. 


w 


APPEHDIX  I 


For  the  reader*  convenience  we  prove  her*  the  result  stated  at  the  beginning  of 
section  2,  concerning  equation  (2.3).  He  assume  that  the  function  v  e  H3,  verifies 
v  •  n  •  0  on  r,  and  assumption  (2.2).  This  last  condition  could  be  replaced  by 
•j|div  vj^  +  l^vl^  <  M.  Let  G  €  Hp,  and  A  be  the  solution  of 

f  -e  AX  +  MX  +  v  •  VX  »  G,  in  0  , 

lee  e 

(5.1)  < 


where  e  is  a  positive  constant.  By  multiplying  both  sides  of  (5.1)  by  AX£>  and  by 
integrating  over  fl,  one  easily  shows  that 


eliX  I*  +  [m  -  (j|div  y\m  +  |7v|„>]  IVA  I*  < 


«  |VGI0  .VXel0  . 


This  estimate,  together  with  (2.2),  gives 


(5.2)  me'o  *  (2/M>  ,7G,0  * 
and  also 

(5.3)  e,AVo  *  ,TC,0  ’ 

Hence,  there  exists  a  subsequence  X£  such  that  Ae  ♦  A  e  H^,  weakly  in  Hg  and 
5  2 

strongly  in  L  ,  moreover,  e  AA^  ♦  0  in  L  .  By  passing  to  the  limit  in  (5.1),  as 
e  ♦  0,  one  proves  that  X  is  a  strong  solution  of  (2.3).  In  particular,  X  verifies 
(2.6).  By  multiplying  both  sides  of  (2.3)  by  X,  and  by  integrating  over  fl,  one  shows 
that 

(5.4)  <M  “  yldiv  v| J  IAI0  <  IGI0  . 

This  gives,  in  particular,  the  uniqueness  of  the  solution  X,  in  H1. 


Sine*  the  linear  map  T  :  ♦  H^,  defined  by  TG  -  X,  ia  continuous  respect  to  the 

2  ""2 
L  no re,  there  exists  a  unique  linear  continuous  map  T  extending  T  to  all  of  L  . 

Clearly,  (5.4)  holds  again.  Furthermore,  X  «  TG  is  a  solution  of  (2.6). 

Remark.  The  result  holds  again  without  assuming  that  v  •  n  «  0  on  r,  and  with 

condition  (2.2)  replaced  by  the  weaker  condition  |div  v|  <  M  (or,  more  generally,  by 

|div  v|s  <  2H)  .  In  that  case,  equation  (2.6)  holds  for  every  text  function  e  e  . 

(we  start  by  proving  the  existence  of  a  solution  X  for  data  belonging  to  the  linear 

space  H  generated  by  an  arbitrarily  fixed  basis  {G  ),  t  ■  1,  2,...,  on  IT)  then,  we 

extend  the  map  G  ♦  X  to  all  of  L^,  by  continuity.) 


1 V V, 


N'.-.’.'-l 

k  *  .•  .  1 


v-’.v 

•  -- 


S  *.**  * 

>  .  v  *  '  » 
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where 

3u  3u  2 

(6.2)  *(u,u)  *  X  I  (-571  +  Tjp)  +  X'tdiv  u) 

i*j  1  3 

Only  for  convenience,  we  will  assume  that  y  >  0,  y  +  v  >  0,  x  >  0,  cy,  X,  X',  and  n  >  0, 
are  constants.  As  for  system  (1.1),  we  Impose  here  the  additional  condition  (1.2). 

The  function  p(p,8)  is  defined,  end  has  Lipschitz  continuous  first  derivatives,  in 
an  l-neighborhood  [m-t,m+l]  x  [n-l,n+l]  Of  (m,n).  By  setting  It  «  pp(m,n), 
y  “  p0(m,n),  one  has  pp(m+o,n+(»)  “  k  -  u^to.o),  pg(m+o,n+o)  »  y  -  i«2(o,a),  where  wi( 
i  ■  1,2,  are  Lipschitz  continuous  (with  norms  <  S)  In  the  l-neighborhood  of  (0,0) » 
moreover  ui^O.O)  -  0.  we  assume  the  physical  condition  k  >  0. 

By  setting 

p  *  m+o,  8  “  n+a  , 

the  system  (6.1)  becomes 

-uAu  -  v7  div  u  +  kVa  +  yVa  ■  p(f  -  (u«V)u]  *  ui^ta.aJVa  +  n)2(a,a)Va  , 
m  div  u  +  u*Vo  +  a  div  u  ■  g  , 

-yAoi  +  yn  div  u  “  ph  -  c^mt-alu*  Va  +  i(u,u)  -  ya  div  u  + 

+  u2(a,a) (a+n)div  u  ,  in  ft  , 

U|r  -  0,  0|r  -  0  , 

the  additional  constraint  is  given  by  (1.3). 


The  linearized  system  is  now 
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